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We investigate a mechanism that generates exact solutions of scalar field cosmologies in a 
unified way. The procedure investigated here permits to recover almost all known solutions, 
and allows one to derive new solutions as well. In particular, we derive and discuss one 
novel solution defined in terms of the Lambert function. The solutions are organised in a 
classification which depends on the choice of a generating function which we have denoted by 



> , 

x{4>) that reflects the underlying thermodynamics of the model. We also analyse and discuss 

00 . 

£NJ ■ the existence of form-invariance dualities between solutions. A general way of defining the 

<N ■ 

• ' latter in an appropriate fashion for scalar fields is put forward. 
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O ■ I. INTRODUCTION 
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During approximately the past three decades there has been considerable interest in finding 
exact solutions to the well know scalar field equations in a flat 4-dimensional Friedman- Robertson- 
Walker (FRW) space-time. This was triggered by the inflationary paradigm [l-4] where a scalar 
field, dubbed the inflaton, plays a central role in producing a brief stage of accelerated expansion. 
The search for exact solutions was to a great extent driven by the need to find a scalar field 
potential which would successfully convey the inflationary prescription, i.e., which would produce 
sufficient inflation, ending with a graceful exit. Another goal was the establishment of a simple 
and concise way of relating the dynamics of the scalar field with the CMB observational data and, 
in particular, to test of the consistency of the slow-roll approximation in the study of generation 



of fluctuations in the large <j) field region 



IB 



The search for exact solutions of models with a self-interacting scalar field has followed diverse 
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strategies. At first the focal point was the obtainment of approximate solutions for the cases 
of potentials assumed to be realistic, because of their particle physics motivation. Examples of 
potentials falling into this class were the spontaneous symmetry breaking double-well potential, 



the potential, the Coleman- Weinberg potential 



Mi. 



Another perspective in the search for 



exact solutions emerged in 1985 by means of which one looked for potentials whose solutions had 
prescribed properties. Lucchin and Matarrese [7] showed the relation between power-law inflation 
and the exponential potential, which was subsequently studied in further detail by Halliwell [§], 
Burd and Barrow [9J, Barrow|lCl|. Ratra and Peebles 11]. Other works along this line explored 



some particular equations of state 
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14j] or arbitrary time dependences of the scale factor [151 ] . 



A third approach has been to explore several recipes for the formal construction of exact solu 



tions. Muslimov [161 ] . Salopek and Bond [13], and J. Lidsey [18( devised a method to obtain exact 
solutions using the cj) scalar field as the independent variable. De Ritis et al Q] explored an 
alternative method based on the Lie symmetries of the equations to derive some exact solutions. 

A skilful combination of the previous methods has lead various authors to derive some classes 
of new solutions 
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26]. 



The quest for exact solutions of single field models proceeds at present, indeed, recently scalar 
field cosmologies have been considered as an explanation of dark energy called quintessence, as 



well as a way to produce the exotic states dubbed phantom matter 
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301 ] . In the latter case one 



envisages the possibility that the kinetic energy of the scalar field be negative. The interest in 



311 ] and to models 
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34]. 



scalar field exact solutions also extends, naturally, to higher dimensional models 
based on modified gravity theories such as, Brans-Dicke or superstring and others 

However two questions have not been fully answered: What is the basic property that allows 
the equations to be integrated? What are the similarities between all the known solutions? 

Here we will be concerned with the question of unifying all the previous solutions in a single, and 
simple framework. The procedure that we have devised, and which is based on a novel generating 
function, unifies all the solutions under a single criterion. Furthermore our method, which goes 
one step beyond earlier attempts to use the <f> field as an independent variable 16l4l8l]. allows us 



to derive new solutions (which we illustrate in subsection IIII B j) . and most importantly provides 
a way to classify all the solutions, and their corresponding behaviours. Indeed not only does it 
enable one to make a complete qualitative analysis of the possible asymptotic behaviour that can 
be expected from single field, scalar field cosmologies, but it also permits us to analyse the form 
invariance dualities that connect diverse solutions and potentials. 

In this paper we work in the framework of General Relativity (GR) and shall consider FRW 
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models with general N spatial dimensions, although many of the applications some of the appli- 
cations will be given for N = 3 to make connection with the literature. First we show how the 
equations of motion for a scalar field ^ in a flat (N + l)-dimensional FRW model can be integrated 
by quadrature. This integration of the equations of motion is not based on any particular tech- 
nique, but is in fact the result of the simple and particular form of the equations of motion of a 
scalar field in FRW space-time. The equations of motion of a self-interacting scalar field contain a 
very special non-linear dissipative term proportional to the square root of the energy of the asso- 
ciated mechanical system. This particular feature is fundamental and underlies the derivation of 
a large of number of exact solutions, and the construction of particular methods to obtain each of 
them that can be found in the literature. In the present work we shall present an unified scheme 
that generates almost every exact solutions for the problem. We illustrate how known-solutions 
are recovered, and, as an example, we shall obtain a solution which to the best of our knowledge 
is novel. 

In the light of the method that we present here, we also analyse the relations that exist between 
apparently disconnected sets of solutions. In this regard we extend Chimento's 35[] and Chimento 



and Lazkoz's [361 ] results, and we discuss a form-invariance symmetry that maps a solution of a 
given set of equations with a given scalar field potential to another solution of different set of 
equations associated with some other scalar field potential. This completes the unification of the 
set of exact solutions of the flat Friedman models with a single scalar field. 

II. EQUATIONS OF MOTION 

Let us consider a (N + l)-dimensional homogeneous and spatially flat spacetime 

N 

ds 2 = -dt 2 + a 2 {t)Y,(dx i f, (1) 

i=i 

where a{t) is the scale factor. Assume that the matter content is a scalar field <fi which is minimally 
coupled. The Einstein equations can be written as 

4> + NH^ + Vtj, = 0, (2) 
±N(N-1)H 2 = ±<P 2 + V(<P), (3) 

where H = a/a, with an overdot representing the derivative with respect to time. Combining these 
two equations one obtain 
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Introducing the new time variable dr = Hdt and the variable 



* = (5) 
we obtain the following planar, autonomous dynamical system 

#' = *, (7) 

where the prime stand for the derivative with respect to r, and is the derivative of V((j>) with 
respect to 0. We study the solutions of © and ([7]) by considering dx/d<p. As 



(iiV(iV -l)-±* 2 ) + ¥ 



x 



we can write 



v 



2x 2xx . 



(8) 



(9) 



N - 1 JV(JV - 1) - x 2 
Equation ([U]) can be approached in two alternative ways. On the one hand, given the potential 
V(cj>) we can, in principle, solve this equation to obtain x ((/)), and thus we can subsequently use (J7J) 
to obtain a solution. But this is only strictly possible for the case of the exponential potential, as 
discussed below. On the other hand, we can instead arbitrarily choose x = x((j>), and in the sequel 
obtain the corresponding scalar field potential, as well as an exact solution through the use of ([7]). 
In this latter case, the integration of equation ([9]) yields 

/ 1 i \ — / x((f))d(f) 

V( ( j>) = A^-N(N-l)--x 2 ( ( j ) ))e N-lJ (10 ) 

where A is an integration constant which determines the amplitude of the potential. The specifi- 
cation of a particular form of x(4>) (satisfying x 1 < N(N — 1)) thus gives the explicit solution 



where one uses equation (J7D, and further utilisation of yields 

1 



H(^) = ±VAe N-1J . ( 12 ) 

Notice that we also have 



, . x(6)d6 

4> 2 = x\4>)e N -U (13) 
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Returning to time t, one gets, 



dt= / jZ- v , (14) 

J ±VAx(</>) exp ( — f x(<fi)d(f>j 

which is the final quadrature. 

By selecting a suitable x = x(<p), and hence H{4>) from equation (fl2j) . one easily integrates the 
equations for the scalar field model in closed form. We obtain the solution explicitly in t cosmic 
time, provided that equation (|14p is invertible. As we shall illustrate in the following section the 
many solutions sparsely found in the literature result from specific, simple choices of x (</>). Our 
result naturally shows how the set of existing solutions can be extended. 

From the definition ([5]) we see that x{4>) is related to the usual barotropic index 7 = (p^+p^) / p<j,, 
where 

p* = Y + v W' (15) 



Y-n<t>), (i6) 



in the following way 



■>'{<■>) Ml — 2 1} 7W> (17) 



and thus choosing x(4>) amounts to a specification of the equation of state. However, the present 
formalism clearly reveals that the choice of the potential V((f>) does not uniquely determine the 
equation of state. Indeed, from equation ([9]), we realise that the specification of a potential is 
associated with families of solutions x(4>) which will differ through integration constants, and 
hence, given (|17|) . will correspond to different thermodynamic regimes. 

If, for instance, one considers x(<p) = A then it follows directly from equation (jlOp that V(<j>) oc 
exp(—2X(p/(N — 1)) (see also Section ITTT|) . and the exact solution takes the form 

N - 1 

<f>(t) = -^—]n(t-t ), (18) 
H(t) = ±^_!)_2_, (19) 

^)^(^4)e-^. (20) 

Notice that the particular case where A = reduces to the case of a cosmological constant, V(4>) = 
AN(N — l)/2 = A, and we recover the de Sitter solution [141 ]. 

But if one assumes instead, from the start, the potential to be exponential (120p . then equation ([9|) 
is integrable, and we derive the most general x((j>) generating function that yields the exponential 
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potential. We obtain though x{(ft) in implicit form 

(/3 - x)~^ a+x \p + xf a ~ x (ax - X) 2aX = r exp (/3 2 a 2 - A 2 ) (21) 



where we have defined a = 2/(N-l) and /3 = ^N(N - 1), and where Tq is an arbitrary integration 
constant. The latter expression cannot be solved to yield an explicit form for x((ft). 

However, we see that we have only two possible choices regarding the limit behaviours of x((ft) 
for the exponential potential for the asymptotic limit (ft —> oo. Consider the change of variables 
X = 1/0 [37], then the planar system ([6|) and ([7]) can be recast in the form 

x' = - (In{N - 1) - \x^j (JL^x - A) , (22) 
X' = ~X 2 x, (23) 

which shows that the infinity manifold x = 0, or equivalently (ft = oo, is invariant. When x = x\ = 
\(N — l)/2 or x 2 = (X2) 2 = N(N — 1), x 1 = also vanishes, and these two solutions correspond 
to the asymptotic, equilibrium solutions. They give completely different asymptotic regimes, both 
dynamical and thermodynamical. A local stability analysis reveals that on the infinity manifold, 
x\ is stable and X2 unstable for A 2 < AN/{N — 1) (otherwise, they will be respectively unstable 
and stable). 

When V(4>) is not exponential we should distinguish the equilibrium solutions of the system 
into those arising at finite values of (ft and those associated with the asymptotic limit (ft — )■ 00. In 
the former case, equation ([7]) yields x = as a necessary condition for a equilibrium point to occur 
at finite values of (ft. Substitution into the equation ([6]) gives 

^.m^in (24) 

and thus the left-hand side vanishes if and only if V $ = 0, i.e., whenever V((ft) has an extremum, 
say (ft = 4>q. A straightforward linear stability analysis reveals that the equilibrium point is stable 
(unstable) if V((fto) is a minimum (respectively, a maximum). 

On the other hand, if the equilibrium point occurs when (ft — > 00, the introduction of x = ^-/4 > i 
as done before, shows that an asymptotic behaviour associated with a non-vanishing constant 
value of x only happens if V{(ft) asymptotes to an positive exponential behaviour, i.e., tv^T^ ~~ ^ ^ 
in equation (|23|) . In this latter case the discussion of the stability previously produced for the 
exponential case applies. 

Still regarding (ft — > 00 case, we also have, of course, the possibility that the equilibrium point be 
characterized by x = when V& — > 0, and hence, V — > A corresponding to the de Sitter behaviour. 
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Looking at equations ([6]) and ([7]) it is easy to see that it is only for the exponential potential 
that these two equations decouple. For an arbitrary potential V = V(4>) this does not occur, and 
thus the condition for decoupling is given by the potential in the form (jlOl) . for a given x = x (</>). 
This leaves us with only one equation to integrate <j)' = x(4>). It seems then natural to assume that 
it should exist a transformation that maps the free-field solution, namely x(4>) = cte, to any other 
solution with a prescribed x (</>). This is the nature of the invariance-form transformation that we 
are going to construct. 

For completeness let us comment on the slow-roll approximation which is a central (if not a 
starting) assumption of almost every scenario in inflationary cosmology. From appropriate combi- 
nations of the Hubble parameter and of its derivatives one can define quantities which take small 
values when the slow roll regime of the scalar field dynamics holds. These are called slow- roll pa- 
rameters, and subsequently the expressions of relevant quantities of the model, such as the spectral 
scalar and tensor indexes, are written as expansions in terms of these parameters in the neighbour- 
hood of the slow-roll regime jf|. Typically, however, only the first few enter into any expressions 
of interest. 

In terms of the generating function x(4>) that we have introduced, the first two slow-roll param- 
eters are 



Apart from a constant of proportionality, e measures the relative contribution of the field's kinetic 
energy to its total energy. The quantity rj, on the other hand, measures the ratio of the field's 
acceleration relative to the friction term acting on it due to the expansion of the universe. The 
slow-roll approximation applies when these quantities are small in comparison to unity, that is, 
when we have both x 2 {4>) <C 1 and §z{(f>) <C 1. This reduces the dynamical system ([6]) and ([7]) to 




(25) 



(26) 




(27) 



4> = x 



(28) 



and this gives 




(29) 



and so ij <C 1 implies that ~ — N, i 



.e., locally V(<j>) ~ exp [-N f x(<£)#] ~ 1 - N f x(<f>)d<f>. 
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III. EXACT SOLUTIONS 



In what follows we list some of the most important exact solutions that can be found in the 
literature, and show how they arise as particular cases of the previous scheme by giving the choice 
of x((f>), which in turn determines the form of the potential V(4>) for each case. This illustrates how 
the present method recovers, and unifies the many solutions that populate the literature. These 
solutions rely on the fact that (I14j) can be inverted, and this together with the fact that in all cases 
x{(ft) can be cast as a logarithmic derivative, x{(ft) = g' {(ft) / g{<ft) where g((ft) is some well-behaved 
function of (ft, are a common feature for all of them. (Since in this section our goal is to establish 
the connection between our procedure and the literature on exact solutions, in what follows we 
restrict to the N = 3 case, and we do not review the details of the solutions. For the latter details 
the reader is kindly referred to the quoted references). 

In the second subsection below, we also derive a new exact solution, and discuss its main 
properties. This enables us to illustrate how the generating method introduced in the present work 
permits to expand the set of known solutions, and is not limited to reproduce them. 



A. List of some of exact solutions 



1. Direct solutions 



1. One of the most important exact solutions is the well known power-law solution which has 
been derived in many forms and through different methods j|. As we have seen in the 
previous section, this solution corresponds to the choice x{cft) = A, where A is a constant, 
and yields (for the N = 3 case) 



V{cft) = A (3 - A 2 /2) e 



\(f> 



(30) 



where A is an arbitrary constant that fixes the height of the potential (in what follows we 
will adopt this notation for the constant factor that defines the amplitudes of the various 
potentials). Note that this solution also permits to integrate the scalar perturbations giving 



a constant scalar spectral index 



2. The Easther solution 



23] corresponds to the choice x(<p) - 
V(<j>) =A(3- <ft 2 /2) e* 2/2 , 



and yields the potential 



(31) 
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where, as explained, A is an arbitrary constant, and subsequently gives 



a{<j>) 



2VA 



(32) 
(33) 



where Ei is the exponential integral function 



381 ] . This solution has the remarkable feature 



that it also yields constant scalar spectral index which is equal to 3. 

The Easther solution is a particular case of the class of solutions characterized by x = \<p. 
In the latter case the solutions are characterized by 



^) = A(3-AV/2) e A <^/ 2 , 



(34) 



and 



ate) = a ^ 1/A , 



i 



Ei 



Ei 



2\ 1 



(35) 
(36) 



where Ei is, once again, the exponential integral function. 



3. The intermediate inflationary solution [12j is given by x(4>) = /V</>, where j3 is a constant, 
and is one of the most famous solutions. The scalar field potential is given in this case by 



Vte) 



16^ 2 ( P 2 
(J3 + 4) 2 V 20 2 



{2Af5) 1 / 2 



-P 



(37) 



and we have 



a(t) = exp [At f 
$ = (2Af3tfy /2 , 

where / is a constant such that < / < 1, and j3 = 4(/~ x — 1) 



(38) 
(39) 



4. In [22] one finds the potentials 



Vite) = A 2 \ 2 [(3 A 2 - 2) cosh 2 te/A) + 2] 



(40) 



and 



V 2 te) = —\ 2 A- 2 <p 2 te 2 + A 2 )te 4 A~ 4 + A - 6 + 2<p 2 ) 



(41) 
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which correspond, respectively, to 

ai (t) = a [sinh(2At)] A2/2 (42) 

<j)i{t) = Aln[tanh(At)] (43) 

Xl {4>) = -2^ 1 tanh(0/A), (44) 

and to 

a 2 {t) = a [sinh(2At)] A2/2 exp[-^ 2 coth 2 (At)/12], (45) 

4> 2 (t) = Acsch(At), (46) 

**) = m 

In the previous epxressions A and A are arbitrary constants. 
5. In [24| one finds a class of solutions with 

<j>(t) = Aexp{-fit n ), (48) 
parametrized by n constant, to which we associate the choice 

mm)^: , (49) 

n 

In Ref. [21(] another family of solutions, also parameterized by constant n, is displayed such 
that 

<j)(t)=A{\nt-B) n . (50) 
this family of solutions is recovered with the choice 

S(<HA) 1 l"eM(<HA) 1 l")<14 

In both cases the general form of the potentials is rather involved, except for some simple 
cases arising form the restriction of n to take some particular values, and thus we refer the 
reader to the cited references for further details about the solutions. 



2. Pair defined solutions 



There are other solutions for which a special method was constructed. The following two 
examples provide a closer look at the procedures that were used. 
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1. In 39j the authors make an assumption which amounts in our prescription to the choice 

[l-F 2 (0)]/3 2 (/3-l) 2 



X[(P) 4-[l-F 2 (0)W-l) 
where (3 is a constant, and thus derived the potential 



V{</>) = exp ^2/3 



IF( 



1 



/ F{<t>) + 17' 

The method provides an exact solution with the further prescription of 

F(<j>) = cosh(A^), 

where A is a constant, then yielding 

V{4>) = A(l + cosh(A0)) T ^ A - 1 . 
Putting (3/X = 2 and choosing the positive sign in ([53|) one gets the solution 



in 



exp(AV^t) + 1 



exp(AV^) - 1 



(52) 



(53) 



(54) 



(55) 



(56) 



2. Reference [40J provides another special method which is equivalent to choosing 



x(</>) = ± 



2g(g) 



(57) 



for which one has 



V(^) = 3H\<p)- g{H{cj>)). 



The solutions of 140(1 were obtained with the choices 



91 (H) = -AH* 



g 2 {H) = ±-VACH-H*H. 



(58) 

(59) 
(60) 

Notice thai ihv authors also derived the solution given in [22] by using a specific choice of 
the function g(H). 

3. In Q a special method was given for a [N + l)-dimensional FRW space-time using certain 
classes of generating functions using 

G(<j>) 



x(<j>) = (N-l)- 



(61) 



fG(<i>y 

where G{(j)) is a function so that G(<j)) = aH{(f>) +L((p), where H(4>) and L((fi) are functions 
specified below, and a is a constant. 



12 



(a) For H-lmeai generating functions and for any N, we recover the solutions of 26] with 
the following choice of x((f>) 



where 



x{(f>) = -(A - l)-^ 1 + a 

L U) e <*4>/(N-D 



(62) 



L(6) = Dae-^^ _J_ e -i^/(JV-i 

^ 3 + 1 



1) 



(63) 



for the cases m = n = 1 and m = n = 2. 
ii. And also L((f>) = ci<p + C2<p 2 when c%, c 2 £ R. 

(b) For the method of the multiplicative generating functions of Ref. 
x((j)) become 

i. G(4>) = oj(fr n and thus 

UJ(p n 



261 ] . the choices of 



x(<P) 



cte + uj4> n+l / (n + 1)' 



ii. G(4>) = uj(j) n exp(— \(p m ) and thus 

, „ e~ x 'f >m m(N- l)(\<j) m y 



x 



(64) 



(65) 



The special cases were considered: n = —3, m = —2 and A = 1/2 for Hq = and 
u = -2/V3; H = 1/2 and ui = 1 + 2/y/3 and with A = 3. 
iii. G(<f>) = ujcp n (l - \4> m Y and thus 

(n + l)(N - 1) (1 - \(f) m ) fl 



x{4>) 



(66) 



The special cases were considered: n = 1, m = 2 and \ = fM = u = 1/2 with 
#o = 1/(24A) and with A = 3. 



4. There are also exact solutions for string motivated models in 
are 

A - 2Be~^ 



4l| . The corresponding choices 



xi{4>) = 2i- 



x 2 {4>) = -2i 



A - Be~^ ' 
A + 2Be~^ 



x 3 (<P) = 2£ 



A + Be-& ' 
A - 3Be~ 2 ^ 
A - Be~ 2 ^ ' 



(67) 
(68) 
(69) 
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where £, A and B are arbitrary constants, which yield the following potentials 



Vi((j>) = A 2 (3-2f)e- 2 ^ + 2AB(Af-3)e~ 3 ^ + B 2 (3-8f)e- 4 ^, 
V 2 {<t>) = A 2 (3 - 2f)e- 2 ^ - 2AB(Af - 3)e- 3 ^ + B 2 (3 - 8f)e~^, 
V 3 (<f>) = A 2 (3-2f)e~ 2 ^ + 6AB(2f-l)e-^ 4 ' + 3B 2 (l-6f)e' 6 ^ 

respectively, and the associated solutions given by: 
1 



H4>) 

H4>) 



A£ 2 
1 

1 

2A? 



-£4>o 



-£</>o 



A - 2Be~^o 
A + 2Be-& 



A \ A + 2Be~^o 

+ 2 n V(l-Ce-^o)(l + Ce-^) 



(70) 
(71) 
(72) 



(73) 
(74) 
(75) 



where C = y/ZB/A. 



B. A new exact solution 



We now derive a new exact solution by exploring the functional forms of (|12p and (|10p uncovered 
by the properties of the Lambert function 42|. Since the solution is novel we revert to the N- 



dimensional case, and later, when appropriate, we shall restrict it to the N = 3 case. 
The Lambert function is defined to be the function satisfying 



W((p)e 



w(<t>) 



(76) 



and it is used in many applications 42|. When cfi is real, for — 1/e < <p < there are two possible 
real values of W((p) [42j. We denote just by W(4>) the branch satisfying —1 < W((f>). Differentiating 
the defining equation (|76|) . and solving for W' , we obtain the following expression for the derivative 
of W: 



W>(<j>) 



1 



(l + W(</>))exp(W(4>)) 

w{4>) 



[i + w{4>)y 



(77) 

(78) 
(79) 



Another useful relation is the integral of W(4>) given by 



W((j))d(j) = (f> W(<f>) + 



1 



W{<j>) 



(80) 
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If x((p) is chosen such that 



-{N -l)W{f{4>)), 



(81) 



then, according to equations ((771) and ((81]) . the pair (V ((f)), H ((/))) that solves the equations of 
motion is given by 

-JV(iV-l) f{4>) {N-lf f\<p) 



V{0) 1 2 W*tf(<l>)) 2 (l + ^(/(0)))2 



from which we derive a first integral 



(82) 
(83) 

(84) 



Another first integral is the Friedman equation ([3]). 

In what follows we consider the N = 3 case. In order to find a solution we take f(4>) = 4> in 
(|8ip and in the following equations. Thus we have x(4>) = —2W'((p) where W(<p) is the Lambert 
function and where A is a constant, In this case one gets for the self-interacting potential V(<fi) the 
expression 

1 W(tf ^ ? 



V ^ A \ 3 2ct>i(i+w((f>)) 2 j w 2 {<t>y 

In Figured] we show a plot of this potential as a function of 



> 




FIG. 1: The self-interacting potential as a function of (j> > — 1/e. 




Using the definition x = <j>/H, it is easy to see that the velocity of the <p field is proportional to 
the condition number of the Lambert function, that is, 

* = ^^ = =^ITWW' (87) 



and using ([80]). one gets, 

± 2VA(t - to) = + . (88) 

These equations cannot be easily inverted, but two different asymptotic solutions can be obtained. 
If (p < e then ([88]) yields 



<j){t) = ±2VI(t - t ) In (±2VA(t - t )J . (89) 

For 4> 3> e one gets 

0(t) = 2V 7 ! * ~ *° . (90) 

W 2TU(±V2/2v^4(t-to)) 
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Because 



W(0)~log0-log(log(0)), 0>1 

the asymptotic form of the potential for large values of (ft is given by 

3(ft 2 1 



V((ft) ~ i4 



_log(0/log0) 2(l + log((/>/log(/>)) 2 
which is thus the form of the potential V for the slow-roll regime. 
The slow-roll parameters 



|, 143] read 

, w 2 (<p) 

'(ft 2 {l + W((ft)) 2 ' 

, w 2 {<p) 

> 2 (1 + W((ft)) 2 



e- 2W ^\2 + W{(j>)) 
1 (1 + W(<ft)) 3 ' 



(91) 



(92) 



(93) 
(94) 



Notice that e < 2 for (ft > and that e — )■ as (ft — > +oo. In Figure [3] we show the values of the 
first and second slow-roll parameters as a function of (ft, notice that both converge to zero when (ft 
goes to infinity. So we see that slow-roll inflation takes place at large values of (ft. At (ft = 0, i.e., 
for small values of (ft, the model has a radiation-like behaviour (since e = 3^/2, and hence 7 = 4/3 
implies e = 2). 

We recall that the condition for inflation to occur is e < 1 and that it ends at e = 1. If we 
assume that the scalar field has a large initial value the inflationary period ends for the value of 
(ft>* that solves the non- linear equation, see ([93]) . 

W'{(ft) = ^. (95) 

It turns out that in this case the value of (ft that corresponds to the end of inflation can be 
determined explicitly. This is one of the most remarkable properties of the Lambert function. In 
order to obtain (ft* consider equation (|78p . Invert, and then multiply both sides by e and notice 
that 



l + W((ft) = W(eV2), 

using (ft = W((ft)e w ^ one gets 



(96) 
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FIG. 3: The first two slow-roll parameters as a function of < 



Let us now study the behaviour of the spectral and tensor indexes in terms of the first order 



slow-roll expansion 



a la 



431 ] . This first order slow-roll expansion is sufficient for our present needs. 



These values of the spectral and tensor indexes are 

n s = 1 - 4(VF'(0)) 2 + AW"{cp) 
n T = -4(lT(0)) 2 . 



(98) 
(99) 



In Figured] we show the plots for both indexes. The smallest value of the scalar field is (j)* in these 
two plots. 



Notice that for large values of 4> we get ns — 1 and ut — 0. Therefore this novel scalar field model 
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FIG. 4: The spectral and tensor indexes as a function of <fi until the end of the inflationary regime. 

based on the choice of x{4>) given by equation (f81~1) . provides another example of a model that yields 
a perfectly scale-invariant Harrison-Zeldovich spectrum for large field inflation. Interestingly also, 
the rolling down of the potential brings the model towards a radiation-like behaviour for small 
values of (p. 

IV. FORM-INVARIANCE MAP 

In this section we address the question: What are the similarities between all the known so- 
lutions? In other words, given that we have shown that the exact solutions are associated with 
different choices of the generating function x((f)), is it possible to go from one solution to another 
solution by means of a transformation x((j)) — > x((p)? The answer to this question is related to 
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the classical technique of obtaining solutions of ordinary differential equations involving a compar- 
ison equation and the construction of a map that performs the transformation from a comparison 
equation to another final equation 441 ], 



The underlying question is whether there is a form-invariance transformation between different 
solutions, that is an invertible map, which we denote by ^, that preserves the form of the equations 

i i n 

of motion. In recent articles, [45j and [35|, this type of symmetry has been discussed. It was shown 
that the Einstein equations of a flat FRW space-time with a perfect fluid 

N(N-l) , , , 
'-H 2 = p, (100) 

p + NH(p + p) = 0, (101) 

where p is the energy density, p the pressure, and H = a/a, admit a non-trivial map, called a 
form-invariance transformation, which transforms the quantities (a,H,p) to (a,H,p), so that the 
latter quantities satisfy the equations 

N(N-l) , , 
>-H 2 = p, (102) 

p + NH(p + p) = 0. (103) 

The ^ map is given by 

* p : (p,H,p) — ► (F(p), ((F(p)/p) 1 / 2 H, -(F{p) + (p/(F(p) f 2 (p + p)dF(p) / dp) , (104) 

where F is an invertible function, that is, 

p^p = F{p) (105) 
H -> H = ({F(p)/p) 1 ' 2 H (106) 
p p = -F(p) + (p/(F(p)) 1 ' 2 (p + P )dF(p)/dp. (107) 

equation (|105p defines the transformation, equation (|106p guarantees the form-invariance of the 
Friedmann equation (|100|) . and finally equation (jlOTj) imposes the form-invariance of the conser- 
vation equation (jlOip . 

If we now consider perfect fluids with a barotropic equation of state p = (7 — l)p, where 7 is 
the barotropic index, the indexes of both fluids are related by 

An application of this form-invariance transformation to the case of a self-interacting scalar field 

rn 

was given in [36] for the case where p oc p. Chimento and Lazkoz have shown that it is possible 
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to transform standard scalar field cosmologies into phantom cosmologies. This phantom duality 
transformation is illustrated by the case of the exponential potential, which as it is well known, is 
associated with p oc p, and constant 7, hence inducing a power-law solution (see equations (|18|) - 
(|20p ). However, no general relation similar to (|104p or f|108|) was given connecting any pair of 
self-interacting scalars fields yielding other types of solutions. 

Here we aim at analysing the general form-invariance dualities that may be established between 
any scalar field solutions. Our formalism is better suited to this purpose than that of reference 36], 



since it relies on the quantities that indeed characterise the scalar field solutions as we have shown 
in the first part of this paper. The quantities considered in (|104p stem from a fluid description 
which does not separate well the roles played by the kinetic and potential energies of the scalar field 
since they are combined in p^. Establishing a form-invariance duality requires both the assumption 
of the transform (1105P and that of an equation of state for the original model, which then yield the 
dual model quantities by means of the Eqs . f|106[) and ([107p . In the case of a scalar field model, 



the method of 36J requires that one be able to integrate back from the assumption of an equation 
of state the corresponding form of the potential, a task which can be cumbersome. In our case the 
form fo the potential is readily available, since it follows directly from the specification of x ((/>). 

We recall that a scalar field can be interpreted as a perfect fluid with the well known cor- 
respondence given by equations (fT5j) and (fl"6l) . and thus it is preferable to use the quantities 
(H, x = 4>/H,c/>) instead of H, p and p. Our aim is to explicitly construct the map 

^ : (H,x,<P) — > (H,x,$), (109) 

such that the equations (|2|) and (|3j) are form-invariant under this map. 

If we let H be transformed into H, the requirement that the scalar field equation, and hence 
the energy density conservation, be satisfied becomes 

% = (5)'<*. (no 

x 2 \H J dH K ' 

In obtaining the latter result the Friedmann constraint equation is also assumed to hold in both 
frames. Another equivalent condition is also derived from the Raychaudhuri equation 

1 dH . . 

= • (HI) 

dH v ' 



In addition, from equation (|T2I) we also get 
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which translates the fact that the form invariance transformation preserves the time variable ([140 . 
and which amounts to be the condition that defines the correspondence between the equations of 
state of the dual models, as characterized by x{4>) and x(4>). 

We establish a form-invariance transformation between any pair of scalar field solutions by 
selecting the corresponding generating functions x((j>) and x((f>), deriving the corresponding H(4>) 
and H(4>) functions in accordance to equation ([120 . and plugging them into equation ([1120 from 
which we derive the relation (p = <f>{<f>). Subsequently, we obtain H = H{H) by using the conditions 
(fTT0[) and (TlTTD . 

Notice that by choosing any two functions x(4>) and x(4>) we are in fact, due to equation (11120 . 
imposing a relation between <p and (p. We will show that by considering x(<ft) = A and x((f>) arbitrary 
it is possible to give an exact solution for this last case in a very simple way. 

A. Proportional Hubble rates 

Let us consider the simple case where the Hubble parameters are proportional, that is, H = cH, 
where c G C\{0} is a constant which can be complex. This case can be immediately tackled by 
resorting to the conditions ([1100 and ([1110 to derive 

= ±Vccj), (113) 



and 



x = ±-^=x. (114) 



so that the form-invariance map is 



^ c : {H, x, 4>) — ► (cH, ±i=x, ±>M), (H5) 



which, given equation ([170 . corresponds to the one given in 35|, [36j. Note that from 

a/a = exp(| , (116) 

one has 

aa~ c = const. (117) 



The previous general equations permit us to review, in a very simple way, the case analysed by 
Chimento and Lazkoz from our viewpoint. In accordance to our results of Section [nil if we choose 
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the pair x(4>) = A and x = X, then from equation (|112p we have 



; e *M<ty = e ^*d0, (118) 



from which we derive 



where 



// () A tfr.A 

— A 



A , <Po (H9) 



7V-1 
6 ° = — 



z/ + In 



H X 2 



(120) 



with f being an arbitrary integration constant. This result subsequently implies that 

We can then distinguish the cases where the constant of proportionality c = (^) 2 is equal to ±1 
from those where it takes some other ratio, and we also distinguish the cases where c is positive 
from those where it is negative. (Naturally the case where c = — 1 requires that one of the A 
parameters be imaginary). Since the time variable is the same in the two solutions which are 
linked by the form-invariance, the cases where |c| ^ 1 correspond to power-law solutions associated 
with different values of the barotropic-index 7, and hence with different equations of state. The 
identity transformation which obviously preserves the equation of state corresponds to c = +1, 
i.e., x = X = X = x. The cases where c is negative are quite interesting since they correspond 
to transformations between standard scalar field barotropic solutions and phantom solutions with 

me case when c = — 1, i.e, A = was discussed in reference 



a negative kinetic energy. T 
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and yields a(t) = a _1 (t) [43, |47]. When c < 0, but c / -1 we have more general transforms 
between standard and phantom power-law solutions such aa~ c = 1. The associated phantom 
form-invariance map 

*_i : (H, x, (f>) — ► (-H, -ix, i(f>). (122) 

gives for any scalar field cosmology its phantom counterpart which is obtained from the map. 

At this point it is appropriate to make two remarks. First, if we do not restrict the generating 
functions x and x to be those leading to power-law behaviours from start, the transform H = cH 
does not necessarily constrain the scalar field model to the case of the exponential potential. 
Second, there is one special application of the proportional Hubble rates transformation this is 
related to form invariance between N + 1 and M + 1 dimensional flat FRW 48J. Indeed, for 



c = ^jp^j ^ one can ma P the solution for (N + l)-dimensional space-time to a 3-dimensional one; 
and for c = -^=1L= from a (N + l)-dimensional space-time to a (M + l)-dimensional one. 
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B. Transforming a power-law solution to any other solution 



Among all solutions described in Section ITU the case of the exponential potential, i.e., x(4>) = A 
stands out as particularly interesting and simple. Letting x(4>) arbitrary it is easy to build the ^ 
map explicitly. Assuming that one chooses the relation between x and x to take the form 



(123) 



</>() 



it follows that if one invert this relation to obtain <j) = f(4>), where A = X/(N — 1), one gets, 



<j> = f{-\- l lnH/H ) 



(124) 



A simple calculation then yields 



-£ = f'(-\- 1 In H/H ), 



and thus 



dH 
dH 



/(—A In H/Hq) 



(125) 



(126) 



Explicitly this map is given by 

9 : (H,x,(j)) — >• | / dH 



/'(-A-MnF/ffo) J - 



(127) 



Notice that by using the solution for the x = A case (|19j) one can construct an asymptotic 
solution for x which can be, in some cases, also a exact solution associated to x. This is given by 

N-l ({N-lf 



m = f 



In 



A ■" V A 2 
In what follows we consider a couple of examples. 



■(t-to) 



(128) 



1. Power-law solutions to the intermediate inflationary solutions 

In the flrs , case we envi S a g e thep„ 88 ibiH, y of mapping the power-law solnticns to the interna 
diate inflationary solution of Ref. [12J], that is the map x = A to x = \/4>. It is straightforward to 
obtain the map 

* : (H,x,<l>) — »• ^^^ 1 - 2 ^,Aexp(-</)),exp(^ . (129) 
Note that /(</>) = exp(</>). 
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Then it follows that the exact solution associated to x = A/0 is given by, using (|128p . 

/ AT l\-(iV-l)/A 

m= -(t-to)-^-^, (i3o) 



(131) 



A 2 

with the potential 

2. Power-law solutions to the new exact solution of subsection 

Let us now turn to the case of x = A and x = —{N — 1)W'(4>). In this case it is straight forward 
to obtain 

and thus 

* : {H,xA) { H ~ HX,{N ~ 2) ^y^)) ■ ( 133 ) 

Then it follows that the exact asymptotic solution associated to x = —{N — 1)W'(4>) is given 
by using (fT28j) . 



with the potential 



w = x (l^_(t^!!W) r « (13 5) 



Expression (|134p extends the solution (j89j) to the N ^ 3 case. 

3. A non-homogeneous transformation 



In 



Parsons and Barrow investigated a transformation which also permits to generate an 
infinite family of solutions for the k = FRW scalar field cosmologies for N = 3. Theirs is a 
particular class of form-invariance characterized by H = a 2 H + /3, and hence cf> = acj), where both 
a and [3 are constants. Applying the equations (jllip and ()112[) of our procedure, we find 

2 = x 2l5Ti a ' ( 136 ) 
aril + p 

so that, after choosing x(cf>) which specifies both H(<j>) and V(<p), given by equations (fTUj) and 
(|12p derive successively H(4>) and x(^). This class of form-invariance transforms illustrates the 
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implication of adding a constant to the original Hubble rate H, that is of a non-homogeneous 
form-invariance transformation. In [49] the authors considered a transform relating old and new 
inflation, and in a subsequent work Barrow, Liddle and Pahud 5(J applied the latter invariance to 
derive a generalization of the intermediate inflation solution. 

If we consider the case of x = A and the map (I136D , which thus reads 



x{4>) = - 



a e -V/(2aO + j3 

we get the corresponding potential, recovering the result of Ref. 



3/3 2 + a 2 I 6/3 



V$) = A 

Also for the = A/ (ft the generating function 

a{4>/a)- x l 2 - 1 



+ 3a 4 e" 



-X<f>/(2a) 



X{<ft) = X 



(/3 + a 2 (0/a)"V2) 



yields the potential 



V{$) = A 



3 - A 



a 2 ((ft/a 



i-A-2 



(/3 + a 2 (0/a)-V2) 



a 



in accordance to Ref. [5 



(137) 



(138) 



(139) 



(140) 



V. SLOW-ROLL, PERTURBATIONS AND FORM-INVARIANCE MAP 

As a joint application of both Sections IIIII and IIVI we discuss the solutions associated with 
the potentials that preserve the slow-roll approximation; this can be done explicitly by imposing 
restrictions on the form-invariance map. 

It is _known that exponential potentials lead to perturbation spectra that are exact power laws 
fl. In 



511 ] a step was given for a systematic classification of types of inflationary potentials that 
yield a constant scalar perturbation indices. The authors obtain the solutions associated with these 
potentials for the Harrison-Zel'dovich case and to general power-laws case both to lower order and 
to next order slow-roll approximation. 

It is possible to think of a infinite hierarchy of expressions for the perturbation spectra and for 
spectral indices. Due to the complexity of the problem, only the first two approximation orders are 
available in general. To obtain the restriction on the form-invariance map it suffices to consider the 
first order of approximation in the slow-roll parameters. The scalar and tensor indices are given 
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by the expressions 



So, using equations ([25]) . we get 



ns — 1 — — 4e — 2ij, 
tit — = — 2e, 

ns — 1 — — 3x 2 + 2a/, 



Since the slow-roll approximation is the dynamical regime where 

e < 1, 
r? < 1, 



(141) 
(142) 



(143) 
(144) 

(145) 
(146) 



it imposes restrictions on x = x{4>). 

Consider x and x, where x corresponds to a slow-roll solution, then using the form-invariance 
map one has 

„2 



X 



2/'(0) 2 

„2 



V 



^-V( 1 + 77j+ 2; ( v 1 -y72j+ x - 



(147) 
(148) 



2 \~'f'J ' ~ \~ j/2 ; 1 ~ y/2 ' 

where / is the function defined by equation f|124|) . This shows that the form-invariance map 
preserves the slow-roll approximation provided that the following condition is satisfied 

1 \ .2 



X' 1 



/ 



12 



< 1. 



(149) 



VI. CONCLUSIONS 



In this work we have presented a unified mechanism that generates exact solutions of scalar field 
cosmologies by quadratures. The procedure investigated here permits to recover almost all known 
exact solutions, and shows how one may derive new solutions. In particular, we have derived one 
novel solution defined in terms of the Lambert function. 

The solutions are organised in a classification which depends on the choice of a generating 
function which we have denoted by x(4>). The choice of the latter reflects the underlying thermo- 
dynamics of the model. Cases in which x{4>) differs only by an additive constant correspond to the 
same potential V((f>). Conversely, this shows that the selection of a potential does not fixes the 
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thermodynamical state of the universe. This is a limitation that must be faced by the efforts of 
reconstructing the potential from observations. 

We have also discussed how one can transform solutions from one class, i.e., characterized 
by a given choice of x(<p) into solutions belonging to other classes. This type of mappings have 



3a]. In the present work we have 



been termed form-invariance transformations in the literature 
extended these transformations to include all sorts of solutions and space dimensions. In particular 
we have generalised Chimento and Lazkoz's results on the duality of standard/phantom solutions 
of power-law models characterized by exponential potentials. We have, for instance, shown how 
one can transform these power-law solution into intermediate inflationary solutions, or into super- 
inflationary solutions either phantom or not. 

We must emphasise that the possibility of establishing a unified procedure to derive exact 
scalar field solutions, and to map different classes of solutions one into another through form- 
invariant transformations ultimately stems from the fact that the field equations are a canonical 
dissipative system in which the dissipative term is proportional to the square root of the energy 
of the scalar field as revealed by the generalised Klein-Gordon equation and by the Friedmann 
constraint equation. 

In forthcoming works we extend the procedure and dualities investigated in the present work to 
exact phantom solutions and to solutions of scalar-tensor gravity theories with a perfect fluid 
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